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We reinvestigate constraints on the eigenvalue density of the Dirac operator in the chiral sym- 
metric phase of 2 flavor QCD at finite temperature, employing the overlap Dirac operator with 
the exact chiral symmetry at finite lattice spacings to avoid possible ultra-violet(UV) divergences. 
Studying multi-point correlation functions in various channels in the thermodynamical limit, we 
obtain stronger constraints than those found in the previous studies that not only the eigenvalue 
density at the origin but also its first and second derivatives vanish in the chiral limit of 2 flavor 
QCD. In addition we show that the axial U(l) anomaly becomes invisible in susceptibilities of 
scalar and pseudo scalar mesons. 
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1. Introduction 

While the classical QCD Lagrangian with Nf massless quarks poses SU(N/)lX SU '(N/)rX 
U(l)yX U(1)a chiral symmetry, U(1)a part is broken explicitly at quantum level by the anomaly 
and then SU(Afy-)z,x SU(iV/)j? part is spontaneously broken to SU(A r /)y in the QCD vacuum. On 
the other hand, it is expected that the SU(./V/)lX SU {Nf)R chiral symmetry is recovered above the 
critical temperature T c , and and the first principle lattice QCD calculations support this expectation. 
It still remains an open question, however, if, how and when the U(1)a part is restored. 

The question whether the U(1)a symmetry is restored or not near T c is of phenomenological 
importance [|]]. Furthermore a connection between the restoration of U(1)a symmetry and the 
gap in the eigenvalue density of Dirac operator near the origin is suggested^]. In principle, the 
fate of U(1)a symmetry can be investigated by numerical lattice QCD simulations. On going 
four simulations with different quark actions, however, have reported different results^], Q ||, [6|]. 
Analytical investigations in the continuum QCD are also inconclusive^, |8|] 

In this report, we address these problems on a lattice, but using an analytic method. We con- 
centrate on Nf = 2 case and employ the overlap Dirac operator [|9], M, which ensures the exact 
chiral symmetry [11 ] through the Ginsparg-Wilson (GW) relation[[l2|] while breaks the U(1)a sym- 



metry by the fermionic measure[[13|]. Using the spectral decomposition of correlation functions 
and assuming the restoration of the non-singlet chiral symmetry, we extract new constraints on the 
Dirac eigenvalue density in addition to the manifest one implied by the Banks-Casher relation^]. 
We also discuss the fate of U(1)a symmetry at T > T c using these constraints. 

In Sec. ||, we explain our setup and assumptions of our analysis. Our main results, constraints 
on the eigenvalue density, are given in Sec. [| The fate of U(1)a symmetry is discussed in Sec. |j. 



For more details of our analysis and results, see Ref. [15] 



2. Setup 

2.1 Spectral decomposition of the overlap fermions 

The quark part of A^-flavor lattice QCD action is given by 

S F = a 4 Y,[YD(A)Y + myF(D(A)) ¥ ](x), F(D) = \-^D, (2.1) 
where \ff = (yi, y^, • • • , y/jvy) r denotes the set of Nf quark fields with the degenerate mass m, a is 



the lattice spacing, and D(A) is the overlap Dirac operator[Q, [Top for a given gauge configuration 
A, which satisfies the GW relation[0] that 



D{A)y 5 + y 5 D(A) = aD{A)Ry 5 D(A) ) (2.2) 

with an arbitrary constant R, and the y 5 hermiticity that D(A)^ = y 5 D(A)y 5 . 

Let us consider eigenvalues and eigenf unctions of D(A), D(A)(j>^ = X% 0„ . The GW relation 
implies that 

A„ A + A,f = aRl^Xt D(A)y 5 tf = XAy 5 <l)A, (2.3) 
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where X A is a complex conjugate of X A , and y^ A is the corresponding eigenf unction. With an 
inner product defined by (f,g) = a 4 Y,xfH x )g( x )> eigenfunctions with complex eigenvalues can 
be ortho-normal as <j>£) = {YsQniYsQm) = ^nm an d (0n >750m) = 0- F° r the rea l eigenvalues, 
X A = (zero modes) and = 2/{Ra) (doubler modes), corresponding eigenfunctions can be 
made chiral, so that we denote the number of left(right)-handed zero modes as N A (N A ) and that of 
doubler modes as n^(n^). 

The propagator of the massive overlap quark for each flavor can be expressed in terms of these 
eigenvalues and eigenfunctions as 



{ImA,?>0} 



fm^+m f m XA+ m 



, f L ti(*)tiW , ^ MMMM! ( 2 4) 

k\ m k\ V( Ra ) 

where f m = 1 —Rma/2, N A +L = N A + N A and = n\ + n A . A measure for a given gauge field 
A can be also written in terms of eigenvalues as 

P m (A)= C - s ™(VAA^ [I (Z 2 m -tftf + m 2 f, (2.5) 

{ImA,f >0} 

where Sym(A) is the gauge part of the action, Ag = 2/{Ra) and Z} n = 1 — m /A|. For even A 7 /, 
P m (A) is positive semi-definite and even function of m. In this talk, we consider the Nf = 2 case. 

It is important to note that all quantities consist of 5a (x,y) and P m (A) are finite at finite vol- 
ume (V < 00), the nonzero quark mass (m 7^ 0) and the finite lattice spacing (a ^ 0). We then 
carefully take the V — > 00 limit and then m — > limit, without worrying about possible ultra-violet 
divergences, until the continuum limit is taken. 

2.2 Chiral Ward-Takahashi identities on the lattice 



With the GW relation, the lattice quark action at m = is invariant under lattice chiral rotation! 1 1 ] 

that 

8a\p{x) = i6T a Y5 [1 —RaD(A)] yr(x), 8 a \j/(x) = ie\j/(x)T a j 5 , (2.6) 

where 8 is an infinitesimal real parameter, and T a for a = 1,2, • • • ,Nj — 1 denotes the generator of 
SU(N f ), and T Q = l NfXNf for U(1) A . 

For the volume-integrals of scalar and pseudo scalar density operators defined by 

S a = a 4 £ [WaF(D(A)) ¥ ] (x), P a = a 4 "£[¥TaY5F(D(A)) ¥ ] (x), (2.7) 

X X 

it is easy to show for Nf = 2 that 

8 b S a = 28 ab Po, 8 b P a = -28 a bS , (for a, b = 1,2,3), (2.8) 
SoS a = 8 a S = 2P a , 8 Q P a = 8 a P = -2S a , (for a = 0, 1,2,3). (2.9) 
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If the chiral symmetry is restored, a product of these operators, & ni ,n 2 ,n 3 ,ns, = P"' S" 2 P$ 3 So 4 with 
a = 1,2,3, satisfies the chiral Ward-Takahashi (WT) identities that 

lim lim T7J (5 a ff n , n2n3tn4 ) m = -5 a0 lim lim =j(JQ&n im m m )m (2.10) 

where 

(0(A)) m = - JmP m (A)0(A), Z = JmP m (A), (2.11) 

J Q = -2*7V 4 L E ^W j 75 ( 1 - jod) ${x) = —2iNfQ(A), (2.12) 

X N=n.k,K V 1 J 

Q{A) = Ng —Ni is the index of the overlap Dirac operator [|l3|] for a given gauge configuration A, 
and / is a minimum integer which makes the V — > °° limit finite in the above equation. Here To 
represents the effect of the U(1)a anomaly. 

2.3 Basic properties and assumptions 

The eigenvalue density for a given gauge configuration A is defined by 

p 4 (A)=limi £ sU-yfifflf), (2.13) 

°° V {lmk£ >()} V / 

which is insensitive to the temperature T, since T is fully controlled by P m (A). Note that p A (A) 
is positive semi-definite for an arbitrary value of X for all A. Although the original eigenvalue 
spectrum at finite V is a sum of delta functions, we expect in the V — > °° limit that p A (A) becomes 
a smooth function. We further assume that p A (A ) can be analytically expanded around A = that 
p A (A) = £~_ p„A n /n! within an arbitrary small convergence radius e. More precisely we here 



assume that configurations which do not have the expansion (2.13) are measure zero in the path 
integral withP m (A). 

In the following analysis, we assume that the vacuum expectation values of m-independent 
function 6 '(A) is analytic function of m 2 in the chiral symmetric phase. Under this assumption, if 
lim„ w o(^ > (A)'°)/m <: = for a non-negative integer k and a positive integer Iq for an m-independent 
and positive semi-definite function GiA), we can write 

>) M = m 2([*/2]+i) j @ A P(m 2 ,A)0(A) l °, (2.14) 

where [r] is the largest integer not larger than r, P(0,A) > at least for some A. In other words, 
the leading m dependence arises from the contribution of configurations which satisfy P(0,A) > 0. 
Using this expression, it is easy to show 

{ff(A) l ) m = m 2 (^ 2 ] +1 ) j $>AP(m 2 ,A)0{A)' = 0(m 2 ^ 2 ^) (2.15) 

for an arbitrary positive integer I, since 0{A) l ° and 0(A) 1 are both positive and therefore share the 
same support in the configuration space. 
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3. Constraints on eigenvalue density 

In this section, we derive constraints on the eigenvalue density of the Dirac operator in the 
SU(2)^x SU(2)« chiral symmetric phase. In this case, the non-singlet WT identities are written as 

Mm lim -j(8a^ ni ,n 2 ,n^ H }m = 0, @n u n 2l n,,n A G ^ (3-1) 

m— >0V— ><x> V 



for o^O, where 



'«l,«2,«3,«4 



n\ +ti2 = odd , ni +«3 = odd , ^tij =N > (3.2) 



is a set of non-singlet and parity odd operators of degree N. Explicitly we have 

~2~^ni,n2>i3>i4 — &n\ — \,n2,njJi4+Y ^2^n\ — 1,713+1,714 ^3 ,«2 + 1 ,«3 — 1 ,«4 ^4^ni+l,«2;"3;"4 — 1 • 

3.1 Constraints at ./V = 1 

At ./V = 1, there is only one operator ^1000 = Pa m &a N , which gives 8 a P a = —2Sq. Using 
the decomposition in Eq. (2A), we have 

lim U-S ) m = lim %{N£ +L ) m +N f (h) m , (3.3) 



where 

for small m. The WT identity now becomes 



lim lim i(-5 )* = lim lim %{N% +L ) m +N f Iim(/i) m = (3.5) 
Since both A^ +L and /i are positive, the above equation gives two constraints that 

lim^« +i ) m = 0(m 2 ), (p A ) m = O(m 2 ). (3.6) 

3.2 Contribution from zero modes at general N 

Before considering N = 2,3,4 cases in detail, we discuss the fate of zero mode contributions 
using WT identities at general N. Let us consider ^i,o,o^v-i S Ua ■ The dominant term in the 
non-singlet chiral WT identity in the large volume limit is given by 

-^) = -(-lfN^/(^ + S\ + 0(V^. (3.7) 

\ / m 

Therefore, the positivity of Ng +L and I\ leads to 

li m ^ N R+Lpm = ( m 2[N/2]+2\ ( g g) 

Since this holds for an arbitrary N, and N^ +L does not explicitly depend on m, we conclude 
that lim (Ng +L ) m /V = at small but non-zero m. 
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3.3 Constraints at N = 2 

At N = 2 there are two conditions from the WT identities, 

X a - n = ^{Sl-P 2 a ) m ^0, ^s = 1_ { s 2 -PX^0, (3.9) 

the former of which has already been examined in the previous subsection. In terms of eigenvalues, 
we have 

Urn*"-* = lim ( -3-Q(A) 2 ) +N f (^+I 2 ) , (3.10) 

v^°° v->°° \ m z V 



in 



where h is defined by 
Using 



-+I 2 = Po- + 2pf + 0(m), (3.11) 
m m 



and (Pq ) m = 0(m 2 ), the condition that hm m ^oX^ 5 = implies 

lun Nf{Q ^P"' =2(pt) m + 0(m 2 ). (3.12) 

3.4 Constraints at N = 3,4 

Five independent WT identities at A = 3 give constraints that 

(Po A ) m = -'^(P$) m + 0(m% lim ^Mk = o(m 2 ). (3.13) 
z v— tn V 

Note that the second condition does not necessary give stronger constraint than (2(A) 2 ),,, = 0(ni 2 V) 
and (Pq ),„ = 0(m 2 ), since it only requires that a set of gauge configuration which satisfy both 
Q(A) 2 = 0(7) and p^ = 0(1) has a weight proportional to m . 

After a little complicated manipulation, WT identities at N = 4 lead toJl~5t] 

(p A ) m = 0(m 4 ), (pf),„ = 0(m 2 ), (pi) m = 0(m 2 ), lim (g( ff )m = 0(m 6 ). (3.14) 

3.5 Final results 

We finally obtain 

lim(p A (A)),„ = (pf) ^ + O(A 4 ). (3.15) 

We think that this condition is the strongest since we know that the Nf = 2 massless free quark 
theory has non-zero (pf )o keeping the exact non-singlet chiral symmetry. Constraints at N = 4k 
(k = 1,2, • • •) in addition lead to (p£ ) m = at a small but non-zero m. 
For the discrete zero modes, we have obtained 

lim \:{Ni +L ) m = 0, lim \-{Q{A) 2 ) m = (3.16) 

at a small but non-zero m, so that these zero-modes give no contribution to susceptibilities in the 
WT identities. 
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4. Discussions: singlet susceptibilities 

For the results in the previous section, the singlet WT identity atN = 2 vanishes as 

= 4 = J im 4^e(A) 2 ), M = o (4.D 

V v-t°°mv 
for small but non-zero m, if the non-singlet chiral symmetry is restored at T > T c . 
For more general cases, the singlet WT identities is written as 

lim lim ^j{5off) m = 2iN f lim Urn ^-AQ(A)0) m (4.2) 

where / is the minimum integer which makes the V — > °° limit finite. After a little algebra, the 
right-hand side of the above equation becomes Jl3|] 



1 _ /QM 2 



Km- t {Q(A)0) m = Km (^-x 0(7°)) =0 (4.3) 
at small but non-zero m, thanks to Eq. ( |3.16 ). We therefore conclude that, for a class of operators 



considered in this report, the U(1)a breaking effect is invisible in the V — > °° limit. This result 
suggests that the chiral phase transition for 2 flavor QCD is likely to be of first order [jl]], contrary 
to the naive expectation that the transition belongs to the 0(4) universality class. 

This work is supported in part by the Grant-in-Aid of MEXT (No. 22540265), the Grant-in- 
Aid for Scientific Research on Innovative Areas (No. 2004: 20105001, 20105003, 23105701,23105710) 
and SPIRE (Strategic Program for Innovative Research). 
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